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We perform a multifractal analysis of the evolution of London’s street network from 1786 to 
2010. First, we show that a single fractal dimension, commonly associated with the morphological 
description of cities, does not suffice to capture the dynamics of the system. Instead, for a proper 
characterization of such a dynamics, the multifractal spectrum needs to be considered. Our analysis 
reveals that London evolves from an inhomogeneous fractal structure, that can be described in 
terms of a multifractal, to a homogeneous one, that converges to monofractality. We argue that 
London’s multifractal to monofracal evolution might be a special outcome of the constraint imposed 
on its growth by a green belt. Through a series of simulations, we show that multifractal objects, 
constructed through diffusion limited aggregation, evolve towards monofractality if their growth is 
constrained by a non-permeable boundary. 


I. INTRODUCTION 

Street networks are ubiquitous worldwide, forming 
possibly the most important and articulated infrastruc¬ 
tural network. These have been considered to be mathe¬ 
matical objects since the eighteenth century [l|, and with 
the discovery of Zipf’s law and the conjecture of Gibrat’s 
law @1 they have gained notable relevance within sta¬ 
tistical physics. In addition, urban patterns share many 
statistical and morphological similarities with biological 
and physical phenomena, suggesting that common organ¬ 
isational principles might be underlying these processes 
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Such patterns can be analysed through fractal geom¬ 
etry. Originally formalized in the 70’s by Mandelbrot 
Q, it wasn’t until two decades later that fractals were 
applied to describe urban morphology [M2 and urban 
growth 00- Until recently, a single fractal dimension 
was employed. It is well recognised now, that a spectrum 
of fractal dimensions needs to be employed to fully char¬ 
acterise systems that present different fractal properties 
at different scales and regions [HI, as is the case of urban 
systems [IM1. These systems are called multifractals. 

Fractal objects are in general described through a sin¬ 
gle measure called the capacity or box-counting dimen¬ 
sion, which measures the amount of space filled by the 
fractal, disregarding local density differences. These dif¬ 
ferences are nevertheless an important aspect to evaluate 
and compare systems according to their local character¬ 
istics, such as measuring performance in cities. In addi¬ 
tion, the box-counting dimension is extremely sensitive 
to external parameters, and has been proven to lead to 
poor results [19,12. 

A measure conveying more information about the in¬ 
ternal components of a fractal structure, directly relating 
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the local density to a measure of proximity between any 
pair of elements of the structure, is the correlation di¬ 
mension 
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Nevertheless, this dimension only gives a general bal¬ 
ance of the distribution of its elements, loosing any infor¬ 
mation on the heterogeneities related to particular zones. 
This is illustrated in Fig.Q] where two very different areas 
are highlighted in the street network of Cardiff. 

The single correlation dimension associated to the 
street network does not give any information on the dif¬ 
ferences between the north and the south layout. The 
latter has more regular and compact structures, while 
the north one tends to be more curly, with a minor num¬ 
ber of intersections. These differences in density convey 
important information about the activities taking place 
in cities. 


In general, high densities of intersections are found 
only in a few places in a city, whereas the majority of 
the network has a small number of intersections. Let us 
denote by <ao the areas with an intersection density larger 
than po, and by aq the areas with an intersection den¬ 
sity pi such that pi < po. The correlation dimension for 
each of the two sets given by <ao and an will be different. 
It is expected that the fractal dimension for the denser 
set (<ao) will be close to two, while that for the set (aq) 
where the points are scattered all over the area, will be 
close to zero [8j|. 

There are many systems that behave in such a way, 
where a single fractal exponent is unable to capture the 
complexity of the fractal structure. This is the basis for 
multifractal analysis , which was introduced by Maldel- 
brot to study turbulent flows [22|, [23|- In general, if an 
object or a process can be described with a single frac¬ 
tal dimension, it is called a monofractal, otherwise it is 
labelled as a multifractal, and an infinite number of di¬ 
mensions could be employed to describe such processes 
[30. Systems that show different local distributions 
for its elements or some of its properties can be commonly 
found in nature, such as the growth probability distri¬ 
bution of a diffusion-limited aggregation process (DLA) 
[2?J, the energy dissipation distribution in a network of 
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FIG. 1: (Color online) The street network of Cardiff, 
UK, along with a zoom of two particular areas. 


resistors [28], the variability in human behaviour [29], 
and the soil’s particle-size distribution [30]. 

In the case of cities it had been extensively argued & 
m, hhh that a single fractal dimension is not enough 
to describe its complex nature. In this paper we perform 
the multifractal analysis of the historical evolution of the 
London street network, by looking at nine digitised maps 
ranging from 1786 to 2010. 

We found that as the city grows, the street network 
progressively fills the available space contained within the 
green belt, and thus creating over time a more homoge¬ 
neous pattern. This process results in the gradual loss of 
multifractality, in the sense that in the most recent time 
layers, the London street network could be described by 
a single fractal exponent. 

We argue that this multifractal to monofractal transi¬ 
tion could be the result of the imposition of the statutory 
green belt around London in 1935 (1938 Act) [33]. leading 
to a condensation phenomena in which the city fills the 
remaining space delimited by the boundaries. We test 
our conclusions by introducing and analysing the growth 
of a diffusion limited aggregation model (DLA hereafter), 
constrained by a non-permeable circular barrier. 


A. Dataset 


In this work we employ a unique dataset, consisting of 
the digitised version of nine historical maps for the area 
contained in the Greater London Authority (GLA), from 
1786 until 2010. This dataset has already been studied in 
several works where further details can be found [l34l-[36| . 

From each of the historical GLA maps, we extract the 
street network corresponding to only the urbanised area 
as defined by the methodology described in [33 • We 
then derive the street intersection map for each of the 
nine points in time, see Fig. [2j Our analysis is then re¬ 
duced to the study of the multifractal properties of street 
intersection point patterns (SIPP) embedded in the two- 
dimensional Euclidean space. 

Though multifractal analysis of point patterns is com¬ 
monly used in many fields (cell colonies, electrochemical 
depositions, galaxy clusters distribution), to our knowl¬ 
edge this approach has never been used before in urban 
studies. We argue that for the statistical analysis, it is 
equivalent to use the street intersections as it is to use the 
street segments. This is the case because the degree dis¬ 
tribution of an urban street network is nearly Poissonian, 
and hence the average degree for the street connectivity 
and its variance are well defined quantities. 

Moreover, the SIPP has been proved to be a quite re¬ 
silient structure over time. Much of the original struc¬ 
ture of the so-called the City of London at 1786 is still 
at place at present times. This situation can be observed 
for almost any major city in the world, which makes this 
approach a very good one to perform a historical analysis 
of cities. 


II. MULTIFRACTAL MEASURES 


A. Review 


In this section we provide a brief review of multifractal 
measure theory. The uninterested reader could skip this 
part and go straight to Sec JIIBi where we summarise its 
applications to street networks. A brief review of basic 
fractal theory is given in the Appendix A. 

Let us recall that when a monofractal object is charac¬ 
terised through a single global exponent, this is computed 
by looking at the distribution density of a specific mea¬ 
sure, and it is implicitly assumed that this distribution 
is uniform [38j. Nevertheless, in real systems there ex¬ 
ist many other intensive quantities of interest, such as 
the temperature, the hardness of a material, the electric 
field, and the density itself among others, whose distribu¬ 
tion is non-uniform, in addition to sometimes presenting 
discontinuities at different scales. 

This led to the definition of a local fractal dimension 
for each of the different regions or subsets of the system 









3 


1786 

1830 

1880 

** 





A ’V. ■% ' 

1900 

m- 

1920 

i ■ 

'f't k' "4 

1940 n 4 

■ *’ -• . V* if i 111 

■' ‘ : tMaBafo 

1965 

y 

1990 

2010 v 

flHSl 

* 

* 


FIG. 2: (Color online) GLA street intersection maps over time. In red we highlight the street intersection maps 
(SIPPs in the text) for the urbanised area as defined by the methodology employed in the text. 


In growing systems, these heterogeneities are observed 
in the different growth probabilities for different regions. 
For example, the DLA has different growth probabilities 
at the tips than at the fjords. These sort of systems have 
been extensively simulated El El through hit probabil¬ 
ities. A probability distribution function can hence be 
constructed from the probabilities pi of growth at loca¬ 
tion i, or from any other normalised measure that sums 
to one in the whole structure. The multifractality of the 
object can hence be described through the distribution 
of the above-mentioned intensive measure. 

Let us denote by /i(cc), the value of a measure fi at 
position x. The amount of this measure within a volume 
V(e,x) in the vicinity of x given by e, is defined by [38j 

H e (x) = / n(y)dy. (1) 

JV(e,x) 

Note that on a typical monofractal structure, this mea¬ 
sure is homogeneous, i.e., 

!*(*)-e~ D . ( 2 ) 

where D is the fractal dimension. 

For multifractal objects the measure p is different at 
different locations, hence the space can be subdivided 
into regions around neighbourhoods sc*, where a distri¬ 
bution for each region can be obtained. If the measure is 


normalised, such a distribution can be constructed using 
Pi(e) as follows 


Pi( e ) — — 


Me (, x i) 

f y(y)dy' 


(3) 


Given the inhomogeneity of the monofractal, and the 
infinitely many singularities of the distribution, for a mul¬ 
tifractal Eq. © becomes m 


Pi(0 


(4) 


where oti corresponds to the strength of the local singu¬ 
larity, and it is referred to as the Lipshitz-Holder exponent 

El. 


This exponent is not unique, as many other boxes or 
subdivisions of the space could have the same oti. The 
number of boxes with the same oti value is given by EH 


y(ai,e) = e 


(5) 


where the function /(a^) is the fractal dimension of the 
set of subdivisions with singularity strength c^. 

To characterise a multifractal system, we therefore 
need to fully specify the measures oti and /(cq). In Fig.[3l 
we give a simple representation on how the different sin¬ 
gularity measures cq of the boxes of size e are related to 
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the fractal dimensions /(c^). Note that if the box has a 
unique measure, the fractal dimension for that box will 
be zero. 
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FIG. 3: (Color online) Schematic representation for the 
distribution of the a — f(a) pair. Each set of has its 
own dimension. In this hypothetical example, the 
dimension A should be the highest, as it covers more 
space (four tiles) than the others. In contrast, the 
dimension C is zero, as it corresponds to a single point. 


Given that we have a distribution function P%{e) for 
the intensity measures in the e-region of a system, we 
can characterise this function via its moments 0 

Z g (e)mJ2PiW- ( 6 ) 

If q = 0, we obtain the number of boxes N(c) of size e 
needed to cover the system Zo(e) = N(e) ~ e~ D . In ad¬ 
dition, via the exponent q we can select different regions 
of the distribution Pi. If q » 0 the regions with the 
largest values of Pi will dominate the measure; similarly, 
if q « 0 the regions with the lowest values of Pi will be 
the ones that will dominate the measure. Given that di 
is the scaling exponent of Pi(e) with respect to e, Z q (e) 
in Eq. © should also obey a scaling relationship of the 
type 

S 9 (e)«£Pi(<0 9 ~(<0- T(,) . (7) 

i 

where r(q) corresponds to the mass exponent and can be 
defined as 


r(q) = qa(q) - f(a(q)). 

(8) 

Eq. ((8]) is usually written in terms of the generalised frac¬ 
tal dimensions D q introduced by 0 as 

r(q) = (q~ 1 )D q . 

(9) 

where 


D q = 1 lim lo ®(^( e ))_ 
q — 1 e^o log(e) 

(10) 


Eq. m allows us to recover three well-known dimen¬ 
sions |44[: for q = 0,Do corresponds to the capacity 
(box-counting) dimension; for q = 1,D\ corresponds to 
the information (Shannon) dimension; and for q = 2,1^2 
corresponds to the correlation dimension. 

If D q is known, the values for a and for f(a) can be 
obtained from Eq. ©, © and m as follows 

a(q) = ±((q-l)D q ), (11) 


q = iL {fiaiq))) • (12) 

Eq. (fl2l) indicates that the maximum of /(a) is given by 
/(a(0)), while all the other regions in the system have 
smaller dimensions. The definition of D q implies D q < D 
where D is the fractal dimension of the substrate where 
our system is embedded, or simply the Euclidean dimen¬ 
sion in this case. 


B. Applications to road networks 

Each one of the multifractal measures has a spatial 
meaning in the context of road networks. That is of par¬ 
ticular relevance in our case, as it allows us to make quan¬ 
titative comparisons between the different SIPP studied. 

Multifractals systems are characterised by a decreas¬ 
ing function of the generalised dimension D q (Eq. (fT0]) h 
while for monofractals D q behaves as a constant func¬ 
tion, and the distribution /(c^) and the measures oti are 
also constants.. 

First, let us focus on the three particular dimensions 
D q , corresponding to q G {0,1,2}. Appendix B outlines 
the methodology step by step to compute these measures. 


1. Capacity (box-counting) dimension: Do 

This has been the preferred quantity identified as the 
fractal dimension, and it gives the probability of finding 
an intersection in an e-box, hence giving a sense of how 
the city occupies the space. 

A high value for Do, for example 1.9, indicates that al¬ 
most all the boxes covering the city contain at least one 
intersection. Nevertheless, this measure is completely in¬ 
dependent of the total mass found in each e-box, which 
means that this dimension tells us nothing about the in¬ 
tersection density distribution of the SIPP. For example, 
two different SIPP could have approximatively the same 
Do while their physical layout might be completely dif¬ 
ferent. 


2. Information dimension: D\ 

For q = 1, we calculate Eq. m when q 1, and 
obtain the following equation: 


D\ = lim 

e —>-0 


-'EjPjlogPj 

- log(e) 


(13) 


where the numerator is Shannon’s entropy, which mea¬ 
sures the amount of information related to the density 
distribution of the SIPP. Eq. (fl3l) is an expression for 
the unevenness of the point density distribution in the 
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e-boxes |44|. Higher values of Di reflect a more uni¬ 
form intersection distribution over the whole city, which 
in turn, could represent either a saturation of the space 
or the fact that the city has a regular structure in terms 
of SIPP. 


3. Correlation dimension: D 2 

When g = 2, we obtain the correlation dimension . 
This is a dimension of the correlations between pairs of 
intersections, II and 1 2, i.e., it is the probability that II 
and 12 lie within the same e-box. 


Do = lim 

e —^0 


log(e) 


(14) 


In this way the correlation dimension gives more ac¬ 
curate results than the box-counting dimension for spa¬ 
tial objects, as it explicitly takes into account the spatial 
structure of the SIPP. A low value for D 2 reflects a spare 
intersection structure, while large values for the same 
quantity represent more compact and ordered structures. 


6. The singularity exponent a(q) 

For different values of g, each region of a SIPP displays 
a different probability Pi of finding an intersection on it. 
As explained in Sec. mi it is possible that different regions 
share the same P^, or at least very similar ones. In such 
a case, the singularity exponent a(q) is a measure that 
captures the variety of intersection densities in a SIPP 
at resolution q. A high value for a(q) means that a high 
variety of densities at that particular resolution is found. 


7. The multifractal spectrum f(a(q)) 

The multifractal spectrum f(a(q)) represents the di¬ 
mension / of the set of regions which display similar a(q) 
values. While the D q represents the dimension obtained 
by examining the distribution of the intersections over 
the whole SIPP, the f(a(q)) is the dimension obtained 
over the different regions which display the same a. The 
plot f(a(q)) vs a(q) gives us a global picture of the whole 
structure of the SIPP in terms of the extension and va¬ 
riety of the different intersection distributions. 


III. RESULTS 


4 . The generalised dimension D q 

For the remnant multifractal dimensions related to 
other q values, we have to keep in mind that they are 
defined via Pi(e) (Eq. (|3])), which is the probability of 
finding an intersection in the i^e-box. 

This differs to measuring the multifractal spectra for 
the DLA |45j, which considers the growth probability on 
the e-box. However, the probability of finding an inter¬ 
section could be related to the growth probability, given 
that the areas in the city with larger Pi have a lower 
growth probability, and vice versa. As stated in Sec. 
ESI q allows us to select different regions according to 
their relative Pi values. In this sense, it can be considered 
a resolution parameter. 

In terms of the SIPP, positive q values correspond to 
multifractal measures for areas with more intersections, 
while negative ones correspond to measures for areas with 
less intersections. 


5 . The mass exponent r 

As stated above, for different values of g, we obtain dif¬ 
ferent probabilities of finding an intersection. The mass 
exponent r(g) reflects this through Eq. ([7]). 

A SIPP with a large r(g) will have more intersections 
at the given resolution q than a SIPP with a smaller r(g). 


A. London street network 

In order to show that London’s street network evolves 
from a multifractal SIPP, to a structure whose multifrac¬ 
tal characteristics are almost lost, we plot and analyse 
the corresponding multifractal measures. To do so, we 
apply equations © to m t° each one of the nine SIPP 
obtained from the historical data set, having considered 
the constraints on the range for the q parameter men¬ 
tioned in Appendix B. 

TABLE I: Main parameters obtained from the 
multifractal analysis of each SIPP. The range for the q 
is not the same for each year due to restrictions 
outlined in the Appendix B 


Year 

q 

Do 

B>i 

d 2 

a(0) 

1786 

[-5.00,14.25] 

1.7959 

1.7697 

1.7467 

1.8287 

1830 

[-7.00,13.25] 

1.7926 

1.7698 

1.7448 

1.8187 

1880 

[-5.75,14.50] 

1.8196 

1.7950 

1.7722 

1.8487 

1900 

[-4.75,19.00] 

1.8434 

1.8252 

1.8109 

1.8658 

1920 

[-4.25,17.75] 

1.8602 

1.8431 

1.8322 

1.8824 

1940 

[-5.25,18.75] 

1.8699 

1.8619 

1.8586 

1.8820 

1965 

[-5.25,16.00] 

1.8850 

1.8803 

1.8780 

1.8927 

1990 

[-4.50,15.00] 

1.8851 

1.8793 

1.8766 

1.8948 

2010 

[-4.50,16.00] 

1.8913 

1.8858 

1.8842 

1.9004 


The actual values for q that have been taken into ac¬ 
count, along with the main dimensions and <a(0), are 
shown in Tab. |H We can observe a consistent asymme¬ 
try between the g + and q~ values, where the range for 
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the positive values is much larger than the one encoun¬ 
tered for the negative values. This means that for each 
SIPP, the number of significant areas with low probabili¬ 
ties of having an intersection in them, which are the ones 
magnified by q~ values, is always smaller than the ones 
with high probabilities. Such a difference is confirmed 
by other measures, and particularly by the multifractal 
spectra a(q) — f(ot(q )), as we show below. 

In Fig. IU panel a, we show D q as a function of q. We 
can see how for the oldest SIPPs the D q values for high 
q are quite low. This relates to the fact that the old 
London street network is scarce with dense areas. The 
increasing of D q for high q during the years shows us 
how these dense areas became more important during 
the SIPP evolution. The evolution of these areas can be 
tracked looking at the distance of D q between q + values 
from one year to another. As the gap between two lines 
is larger, the morphological difference between structures 
is also larger. As an example, for q = 5, the SIPP dis¬ 
tribution between 1880 and 1900 is less similar than the 
one between 1900 and 1920; the differences for all q val¬ 
ues between the last three years is almost zero, i.e., it is 
not trivial to distinguish one from the other in terms of 
the distribution of its intersections. This is an evidence 
of how the London’s SIPP evolved in time toward a con¬ 
densation, as we better explain below. The inset plot at 
Fig. 2J panel a, shows the evolution of Do, D\ and D<± 
during time. The clear difference between these values 
(Do > D i > D 2 ) from 1786 to 1920 is a clear signature 
of the system’s multifractality during this period. From 
1940, the nominal distance between D\ and D 2 begins 
to decrease, until 2010 when is practically zero ( 0 . 001 ), 
while the distance from Do and D\ at this same year, is 
only 0.005. This strongly suggests a decay in the SIPP 
multifractal characteristics over the years, almost to the 
point to be able to describe the last three SIPPs with a 
single generalised dimension, ie, D 2 . 

Such a phenomenology is related to a homogenization 
process for the SIPPs, possibly related to a condensation 
phenomena, as we show below. 

The mass exponent r(q) as a function of q has the 
expected monotonic growing behaviour, as the SIPP 
number of intersections increases year after year. Until 
q m 3, r(q) is practically the same for all years. As the 
q values increase, the areas with more intersections are 
the ones enhanced, and it is there when the differences 
between r(q) values emerge. From Eq. (j 8 j), we know 
that t(0) = —D q , t( 1) = 0 and r(2) = D q . Then, as 
it has been argued, different values for D q imply multi¬ 
fractality. As the non-linear behaviour for r(q) would, 
in turn, imply multifractality, we test the linearity for 
the r(q) for each year, fitting it against a linear function 
and looking at the norm of the residuals of such fit. 
If the norm is equal to zero, then the function can be 
characterised as linear. For the first three years, we have 
higher norms (all > 1 ), but as we move in time, this 
norm is approaching zero, a situation that implies lost 
of multifractaliy during the growth of the London SIPP. 


In panel b of Fig. 0J we show the multifractal spectrum 
f(a(q)) as a function of q, for the SIPPs. The function 
f(a(q)) conveys information with respect to the distribu¬ 
tion of the intersections, in the sense that it relates the 
different areas enhanced by q with a dimensional quan¬ 
tity. 

We observe that for q~, the curve follows the same 
behaviour for all the years, indicating that the number 
of areas with low probability of having an intersection 
remains more or less constant during the years. 

On the other hand, we observe significant structural 
differences in the SIPP between the different years for 

q + ■ 

In particular, we observe a change in the concavity of 
the curve once we reach 1920, and for the most recent 3 
points in time, the curves almost overlap at q ~ 5. 

Let us denote by k the mean curvature of f(a(q)) 
in a closed neighbourhood of this point. f(ot(q)) has 
a maximum at q = 0 , and for a monofractal structure 
f(a(q < 0)) = f(a( 0)) = f(a(q > 0)), hence k = 0. For 
a multifractal structure, we expect to find k 7 ^ 0. The 
plot tells us how the value of hi decreases with the evolu¬ 
tion of the network, reaching a stable low point for the 
past 50 years. 

In panel c of Fig. [H we show the behaviour of the sin¬ 
gularity exponent a(q) as a function of q. It represents 
the structure of the SIPPs in terms of the dominant prob¬ 
abilities at resolution q. Once again we observe that for 
q~ the structure of the SIPP is very similar for all years, 
while for q + major differences between the years emerge. 
For the earlier years, there are only a few areas with ma¬ 
jor intersection densities, and so the values for a(q) are 
low for q + . On the other hand, for the latest years, this 
situation is inverted. 

From q = 0, the values for a(q) for the three most 
recent SIPPs are very similar, i.e. the distribution of 
probabilities is the same for all areas (less/more inter¬ 
sections), while for the rest of the SIPPs, the differences 
between these values are greater, particularly for the first 
two ones. 

Nevertheless, from q = 7, each of the values for the 
SIPPs tend to converge to a constant value, indicating 
that from this point the different spectra remain constant 
and hence the resolution limit for q has been reached. 

In panel d of Fig.0J we show the multifractal spectrum 
/(<a(g)) as a function of a(q). 

As we already mentioned, each point of the plot shows 
the fractal dimension f(a(q)) of the areas with the same 
singularity strength a, at resolution q. In this sense, it 
is important to note that f(a) does not represent the di¬ 
mension of continuous regions, which tends to be a com¬ 
mon misunderstanding. We observe that these curves 
have a maximum for different a(q) values. These maxi¬ 
mum values correspond to Do , at <a( 0 ). 

It might then become clear that the a(q) values to 
the left of the maximum are the ones associated with 
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q > 0, while the ones to right are the ones associated 
with q < 0. Keeping in mind that D q is a dimension 
for the whole structure, while f(a{q)) is the dimension 
of a subset of such structures, the clear asymmetry of 
the curves reflects at some extent the asymmetry in the 
distribution between regions with large/small number of 
intersections. Then the differences between the spectra 
for areas with high density of intersections are greater 
than the differences we find for the areas with small in¬ 
tersection densities. 

The fact that the dimension f(a(q)) is greater for the 
zones with a values associated with q < 0, which are the 
areas with less number of intersections, is not a surprise. 
For example, for 1786 the f(a(q)) dimension for the areas 
with larger number of intersections is very low; while the 
f(a(q)) dimension for the areas with less intersections is 
larger. 

If the multifractal spectrum collapses to a single point, 
it means that the dimension f(a(q)) becomes the same 
for all ct(q), i.e. either all the regions in the structure 
share the same strength exponent or all the different 
strength exponents have the same dimension at all reso¬ 
lutions. In either case, this means that the structure is a 
monofractal. 

For the case of the SIPPs, it is clear that none of them 
collapse in such a way, nevertheless we can observe a 
transition from a wide curve in 1876 to a narrow one for 
the last two SIPPs. This indicates a progression from 
multifractal to mo no fractal characteristics. 

The symmetry of the curve a(q) vs f(a(q)) gives us 
information about the relationship between regions with 
different distributions. In panel e of Fig. 0J we show a 
symmetry measure for the multifractal spectrum. This 
is done by plotting the points {|a(0) — a(q^~)\, |<a(0) — 
a(q +)|} for each year, where qi and qj stands for the 
last/first valid q value, and the 1:1 correspondence line. 
The distance from these points to this line is a visual tool 
to understand the asymmetry of the relative spectrum. 

This symmetric feature represents the balance between 
areas with more intersections and areas with less inter¬ 
sections. 

The actual distance for each point is shown in the 
inset of this same plate. It is clear that this distance 
increases as we move forward in time: for 1786 the 
point lies at a distance of 0.02, while for the last three 
years the distance ^ 0.17. This loss of symmetry, 
due to the increasing number of intersections, reflects 
the undergoing condensation phenomena in our study 
area, and once more, supports our hypothesis on the 
evolution of London from a multifractal to a monofractal. 


B. Model: Constrained DLA 

The multifractal measures show that London’s street 
network grows while undergoing a process of morpho¬ 


logical re-structuring that is manifested in fundamental 
differences in the distribution of the structures between 
different time periods. The behaviour of the network, 
with respect to the way it fills the space, follows an ex¬ 
pected trend until 1920, and from this year an anomalous 
evolution seems to take place, with D\ and D 2 starting 
to collapse to a single value. We argue that such a change 
comes as the result of two factors: the extensive diffusion 
of development from the core of the city to the periphery 
with massive suburbanisation and, after the mid-1950’s, 
the effect of the actual implementation of a green belt 
around London, in order to contain urban sprawl. 

In previous work we showed that the city undergoes 
a condensation process in its evolution as the result of 
the green belt j3J , here we show its effects on the inner 
structure of the city for the first time. 

In order to test our hypothesis of the evolution of Lon¬ 
don from a multifractal to a monofractal structure due 
to the presence of the green belt, we introduce a mul¬ 
tifractal diffusion model in a constrained plane, and we 
explore the effects of a restricted barrier on a growing 
multifractal structure. Our model is based on the clas¬ 
sic DLA algorithm [27] modified in a way that mimics a 
green belt barrier set at a distance d from a seed point 
at (0,0). 

Given that the DLA structures are formed through a 
stochastic process, we perform 100 runs to derive some 
average measures, and since the historical data set used 
in this investigation consist of nine points in time, we 
construct a nine point time series T for our model. Each 
point Ti G T is selected at one particular iteration, rep¬ 
resenting different stages of growth. 

We define the barrier as a non-permeable circular perime¬ 
ter c, with radius r = 100 and centre at (0,0). The model 
follows exactly the same rules as in the DLA, until the 
wandering particles enter in contact with the green belt. 
At this stage, the particles are not allowed to cross the 
barrier and are forced to find a new position in the avail¬ 
able space inside the barrier. 

The effect of the barrier on the DLA model can be ob¬ 
served from the sixth stage of growth onwards, see Fig. [5j 

We then apply equations (j6j) to m to each aggregate 
in order to investigate the effects of the growth restric¬ 
tion on the model with respect to the inherent multifrac¬ 
tal characteristics of the DLA (see Fig. El )• We observe 
that all the measures related to multifractality start to 
collapse toward monofractality after the sixth point. In 
the last two stages in Fig. [5j the model is practically 
a filled 2-D circular plate, with a correlation dimension 
D 2 = 2.0. 


IV. CONCLUSIONS 

Several studies have already showed that cities are mul¬ 
tifractal objects, instead of fractals structures, as it was 
originally suggested. Here we show that for London, this 
is only true up to a certain point in time, when the city 
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FIG. 4: (Color online) Multifractal Analysis for all the years studied. The differences in value between these 
measures follow a pattern that imply that a single dimension is enough to characterise the SIPP, or, in other words, 
the structure is undergoing a transition from a multifractal to a monofractal structure. 


starts a condensation process forcing the new street in¬ 
tersections to be allocated in the available space inside 
a confined area, i.e., not in a sprawling fashion. For the 
three most recent maps examined in this analysis (1965, 
1990, 2010) all our measures are very similar, reinforc¬ 
ing the idea that a condensation process is in place and 
that it is reaching a stable phase. This implies that oti 


in Eq. (ED is equal V g, i.e. the probability of finding an 
intersection at any random area in the SIPP is the same. 

We observed that such a multifractal to monofractal 
transition is not particularly smooth, in the sense that 
the differences in values get dramatically shortened from 
1940 and 1965, which is the period when the green 
belt was established de facto. After 30 years of the 
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FIG. 5: Evolution of the constrained DLA model 


introduction of the green belt, the structure lost part 
of its multifractal signature, while after 60 years we 
can practically characterise the city as a monofractal 
structure. 

In order to investigate the impact of such a natural bar¬ 
rier on the internal morphology of the city, we performed 
a simulated experiment. We imposed a non-permeable 
barrier on a well-known multifractal diffusion process 
(DLA) and measured the evolution of its multifractal 
properties as the structure grew and began to approach 
the barrier. We obtained a clear transition from multi¬ 
fractal to monofractal behaviour in the two final stages of 
the growth process of our structure. This is in very good 
correspondence with the observed behaviour for London. 

The results presented here are hence another step to¬ 


wards the understanding of the implications of restricted 
urban growth policies in the urban morphology of mod¬ 
ern cities. 


Appendix A: Fractals review 

Mono fractal objects are characterised by invariance 
under a change of length scale [8, TjJ. When the change 
is performed in an isotropic way, the fractal is called self¬ 
similar. If on the other hand, the length scale is changed 
in different directions through different factors, the frac¬ 
tal is called self-affine. 

The scaling exponent of the power-law relationship be¬ 
tween the size of an object, e.g. its mass M, and the ob¬ 
servable being measured, e.g. the length L, is in general 
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FIG. 6: (Color online) a) q — D q plot. Until point 6, 
when the structure begins to interact with the green 
belt, all the points for both models depicts the same 
similar behaviour. From point 7 the effect of the barrier 
over the structure is clear, b) a(q) vs f(a(q)) plot. The 
spectrum at points 8 and 9 are almost completely 
collapsed to a single point, which is an indicator of a 
monofractal structure. 


identified as the fractal dimension D [47] 

M~L d . (Al) 

In this sense, as the fractal grows, its density p measured 
in the Euclidean space d decreases: 

a = d — D. (A2) 

It is important to note, that this fractal dimension is not 
always unique. Different ways to characterise the system, 
and to measure scale-invariance with respect to different 
properties might lead to different exponents. For exam¬ 
ple, if a metric is defined in the space, one can find the 
scaling relationship between the number of units needed 
to cover the space as the fractal object is re-sized. If those 
units are balls, the exponent in this case corresponds to 
the Hausdorff-Besicovitch dimension Dh 0, Ell • This 
reduces to a problem of optimization of the number of 
balls needed to cover the space. In the case of simple 
self-affine fractals, the different fractal dimensions con¬ 
verge to a single value, and the fractal can be fully char¬ 
acterised by a single dimension D = Dh • Nevertheless, 
this is seldom the case for fractals encountered in nature. 


a. Random Fractals 

Fractals are classified in terms of deterministic and ran¬ 
dom fractals [50]. Deterministic fractals correspond to 
mathematically constructed objects whose scale invari¬ 
ance holds for all scales. Examples of these are the well- 
known Sierpinski gasket and the Koch curve. In experi¬ 
ments and in nature only random fractals are observed. 
These are finite objects that can be characterised as frac¬ 
tals within a specific regime of length scales. We denote 
by a the smallest linear size that can be measured, and 
by L m the largest size, that might correspond to the size 
of the box where the finite fractal can be embedded in a 
Euclidean d-dimensional space |47]. For random fractals 
Eq. (jA 2 j) becomes [Meakinl990]: 

p(X) ~ A-“; A = —. (A3) 

a 

In general the linear size L should be replaced by L/a. 
At this point it is very important to note, that random 
fractals exhibit self-affinity or self-similarity symmetries 
at a statistical level only. This means that the correct 
way to measure the scaling exponent is through a large 
number of samples. The density ceases to be global, and 
one needs to take into account the average over all origins, 
and for non-isotropic fractals, over all orientations and 
ensemble realisations. The density is hence replaced by 
a density-density correlation function 

°( r ) = y^p( r + r ')p( r ')- ( A4 ) 

r' 

that can be seen as the probability of finding a particle 
at r + r' if there is already one at r', and where the vol¬ 
ume V can be considered in terms of the total number of 
particles. For isotropic fractals, the density-density cor¬ 
relation function corresponds to the density distribution 
around r, since c(r ) = c(r'). This is scale-invariant 

c(rjr) ~ r]~ a c(r). (A5) 

where 77 denotes the re-scaling factor, and can also be 

expressed as 

c(r)^r -a . (A 6 ) 

For growing systems whose mass can be characterised by 
the number of points (or particles), the fractal dimension 
can be obtained by looking at the scaling relationship 
between the number of points M(R) within a sphere of 
radius R [51] 

r R 

M(R) ~ / c(r)d d r ~ R d ~ a ; D = d - a. (A7) 

Jo 

Note that this equation only holds at the asymptotic 
limit, i.e. L m and R a. In addition, this radius can 
be seen as the radius of gyration , and the system can 
be considered as a growing system as presented above, 
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or as a system that contains several fractals of different 
sizes that are the outcome of the same process. For this 
latter case, the average needs to be considered, and in 
order to reduce statistical fluctuations, a large number of 
samples is required. In a different paper we will investi¬ 
gate this scenario, considering cities in a country as the 
different realisations of the same process. The observed 
fluctuations are the result of the statistical nature of the 
self-similarity of these systems. 

Obtaining a good estimate for the fractal dimension is 
not an easy task jUj. For small scales of L, corrections 
to the above given equations need to be introduced. In 
addition, further inaccuracies arise from the limitations 
of the sample size and the orders of magnitude involved 
in the distribution of fractals from which D is estimated. 

One solution is to verify that all such fractals involved 
in the computation of D are self-similar. This is done 
by collapsing all the curves for the different fractals to a 
single one using the equation for self-similarity Eq. (1A5[) 
and by taking the characteristic length of the system as 
the scaling factor 

(A 8 ) 

R 

leading to 


c(r/R) ~ (4) c(r), 

(A9) 

c(r) - R~ a f{r/R), 

(A 10 ) 

~ M^ D ~ d)/D f(r/M 1/D ). 

(All) 


b. Hurst exponent 

The fluctuations of a system around a specific pattern 
can also be characterised in random fractals that are self- 
affine through the Hurst exponent H [H, 45|. Self-affine 
fractals present symmetries that are not isotropic, lead¬ 
ing to different scaling factors in different directions. In 
a 2 -d space, the coordinates in the x and y axes might be 
re-scaled via r\ x and y y , such that: rj y = rj x - The above- 
mentioned fluctuations can also be observed in time se¬ 
ries, and H will give an indication of the randomness of 
the data, or in single valued functions. If time for exam¬ 
ple is re-scaled by b (for b > 0), the whole function F(t) 
needs to be re-scaled by a different factor 

F(t) ~ b~ H F(bt). (A 12 ) 

An important example is the Brownian motion of a par¬ 
ticle [Meakinl990]. Let us denote by B{t) the distance 
covered by a particle in time t. In this case, if t bt, 
then the process B{bt) is invariant if and only if a change 
in the distance scale of b 1 / 2 is introduced simultaneously: 
^^Bit) ~ B{bt). Hence H = 1/2 in this case, confirm¬ 
ing self-affine symmetry. 

Other important examples are the financial time series. 
There is a lot of controversy with respect to the scaling 


properties of such time series p53l455| and there are views 
that such processes should be considered as multi-scaling 
processes, see [56[ for a review. 

Appendix B: Practical calculation of the multifractal 
measures 

The formalism introduced in SecJlT] presents several 
practical disadvantages. The two most important ones 
are: 1) the smoothing of the curve D q , and 2) the nu¬ 
merical calculations of the derivative in m and m- 
To overcome these difficulties, we apply the method de¬ 
scribed in [E3 to obtain the parameters D q , r(g), a(q) 
and / (a(q)). 

We proceed as follows: 

1 . Calculate the number of non-empty boxes of size 
e,N(e), necessary to cover the SIPP. 

2 . Calculate the number of intersections at the box 
i of size e,/i(i,e), which is the discrete version of 
Eq. (1). 

3. Calculate P(i, e) given in Eq. (4) 

4. Calculate the partition function of Eq. ( 6 ) : 


Z(q,i,e) ~ 

(Bl) 

5. Apply the normalised measure //,(</) 


mHH) = r 7( . V 

z(q,ht) 

(B2) 

6. Obtain by linear regression, a(q, e), /(g, <a, 
r(q,e) : 

e) and 

a(q,e) = fi(q, i, e) e)). 

(B3) 

f(q, a, e) = i, e) In(M<?, h e))). 

(B4) 

T(q,e) = X?Jl ) P(i,ey~ 1 . 

(B5) 

7. Finally calculate D q according to 


C( ? ,«) = r(? ' e 1 ) . 

(B6) 


In practice, q does not take values in the entire range 
(— 00 , 00 ), due, first, to the obvious computational lim¬ 
itations, second, to the inherent statistical errors asso¬ 
ciated with the linear regressions and third, to the spe¬ 
cific dimensionality constraints of our SIPP. The topolog¬ 
ical representation of these networks is a two-dimensional 
point set and Eq. m implies that /(<a(g)) < D, where 
D is the dimension of the substrate on which the measure 
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is distributed (in our case D = 2). But Eq. (flQl) can in¬ 
deed generate values D q > D and f(cx(q)) > D for q < 0. 
In this research, initially we selected q in the discrete in¬ 
terval [—20, 20] with a 0.25 step separation between each 
values, and the actual valid range is selected in two steps: 

a: All the q values where D q and f(a(q)) are greater than 
2 are dismissed for the further analysis and, 

b: For the values obtained through a regression, we se¬ 
lected the ones such that R 2 > 0.9. 
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